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AA 


Capsule Entry 


January 15, 2006 


sie} [e[=) 


Stardust 


(Figure from http://appel.nasa.gov/) (Figure from http://en.wikipedia.org/wiki/Impact_event) 


¢ Can some of the modern computational analysis tools used in design of 
heatshields of capsules be used (repurposed?) for simulation of meteoroid 
entries? 
¢ Can we build or develop, across various classes of meteoroids, models for: 
¢ Material thermal response 
¢ Material structural response, including fragmentation 
¢ Energy deposition along meteor trajectory in the atmosphere, /.e., light curves 
¢ How much would the results of these models differ from, and improve upon, 
those obtained from the equations of meteor physics? 
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Entry Capsules vs Meteors eS 
Meteoroid (Asteroid) 


A - 
Ames Research Center 


Shape Regular and smooth geometry Irregular and rough geometry 
Material | Manufactured ablative material | Depends on asteroid class (S, M, X) 
eienite | Minimal voids and cracks; | Could have voids & could be fractured; 
Known structural properties Structural properties of meteorites 
Shape Change | Negligible recession/mass loss | Recession dominated 
Fragmentation Not an option! S-class most likely to fragment 
Recovered panei “i 
Stardust i dene 
Meteorite 


Capsule 


(http:/Awww.nbcnews.com) a 
¥ 


Capsule (Earth entry) Meteoroid (Asteroid) 


Validation Data Sensors, imaging, capsule recovery | Light curves, infrasound, falls 


Crounculectin Ballistic ranges, arc jets, wind Probably same facilities as capsules, 
g tunnels, laser facilities but with limitations 


Although meteor physics has some things in common with capsule entry physics, 


the approaches to the problem are different — reconstruction vs prediction 
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Focus Areas of Modeling 


Ames Research Center 


¢ Flight Mechanics 
¢ Ballistic vs Lifting entries 
¢ 3DOoF vs 6DoF 
¢ See poster by Allen in this workshop 


Aero/Aerothermodynamics 
¢ Turbulent convective heating 


¢ Thermal radiation heating 
¢ This presentation 


e Material Thermal Response 
¢ Ablation & Recession 
¢ Melt vs Vaporization 
¢ See poster by Stern et al. in this workshop 


¢ Structural Response 
¢ Pressure & shear loads 
¢« Dynamic loads 
¢ See poster by Agrawal et a/. in this workshop 
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A - 
Ames Reseorch Center 


Mass Loss Model — Single Body eS 


p, Ambient density QO Heat of ablation 
u,, Meteor velocity m,, Meteor mass 
C,, Heat transfer efficiency t Time 


A,, Area (x-section or wetted?) 


LHMEL Test of H chondrite 
: i : (Tamdakht) @ 5 kW/cm2 
¢ C,, is efficiency of conversion of freestream energy 
into heating of meteor surface 


¢ Heat of ablation, Q, is a big source of uncertainty 
e Need to understand energetics of melt vs vaporization 


¢ Exploratory test on meteoritic materials performed at 
LHMEL* 


¢ Surface irradiation with CO, laser 


“Laser Hardened Materials Evaluation Laboratory 


(see poster by Stern et al. for details about estimating Q) 
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‘What can flow computations provide? Guy 


Ames Research Center 


e Flow computations for a (hemi)spherical geometry can provide C,, 


Flow computation Radiation computation 


(DPLR code) (NEQAIR code) 
2 20R-q(O)sin@ dé . a 
oe Ca ee i WOV=Gi ae OP Gira ne @) 
3 
—o uA 
2 Pa, mm 
Radiative 2 Radiative Tauber-Wakefield 
Radiative 0) = adiabatic -_ Adiabatic ‘gs . 
eae eee | 1 3 AT 1 radiation cooling 
(1+3.41) 5 Pall correlation 


e Flow computations can also provide estimate of energy radiated in a 
specific wavelength interval, e.g. V-band, from the shock layer 
e Luminosity can be converted to a “magnitude” 
¢ Construct a synthetic light curve for direct comparison to observations 


Pos = 22m, sind| f° 12", (A:0)d2)0 
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Flight Space Concept" 
(Uses US 1976 Standard Atmosphere) 


2 Chelyabinsk data from Borovicka et al., Nature, 2013 ‘ 
S > Scaling Laws 
o Stardust" P.tag/bar 
cs * | 7 
1 
w= Conv 3 
© “x ——./O Uu 
~ D stag Pa m 
gq Oo JR 
oO m 
Ss oO 
5 SI Chelyabinsk oe x R p,u.. 
= (Main Fragment) 
© 
~N 
Tropopause © j ij = 
popau Flight Space Hemispheres of radii (R,,) — 1, 3, 10, 
f-) 20, 30, & 100 m 
0 5 10 15 20 25° Spheres of radii— 1, 3, 10, 20, & 30 m 


Velocity/km.s- 
e Flight Space is a way to delink trajectory from environments; X-33, Orion, ... 
e Flight Space is specific to a given shape — (hemi)sphere for now 
¢ Compute aerothermal environments at nodes 


e Heat transfer efficiency and energy deposition are functions of R,,, 0,(Z), u,, 
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2 a 
A High-Fidelity Computations vasal 


A _ 
Ames Reseorch Center 


The meteor body does not ablate and does not cool by re-radiation (cold wall) 
e Allows application of physically meaningful surface boundary conditions 
e Provides the upper bound on heating (convective and radiative) 
¢No blockage by vapor phase of meteoritic material!! 


e Flow computations: 
e Axisymmetric Navier-Stokes calculations for body in a fixed frame of reference 
¢ Turbulent flow of 11-species air (N5, O., NO, N.*, O.*, NOt, N, O, Nt, Ot, & e-) 
Does not account for any surface roughness 
eGas phase rate chemistry, but thermal equilibrium 


e Radiation computations: 
e Decoupled from flow computations (adiabatic inviscid shock layer assumption) 


e Line-by-line simulations with temperatures & number densities from flow solutions 
Includes discrete transitions (atomic lines) and continua (bound-free & free-free) 
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Convection vs Radiation 


Sample surface heat flux distribution 
(u,, = 20 km/s, p,,,. = 30 bar, R,, = 15 m) 


—Convective Radiative Total 


— 
oO 
WO 
— 
© 
N 
a 
0 10 20 30 40 50 60 70 80 90 


Angle from stag. point/deg 


Ames Reseorch Center 


Heat flux/kW.cm2 
400 


¢ Surface heating completely dominated by shock-layer radiation 


e True across all velocities and hemisphere diameters, except for small (1 m diameter) 
hemispheres at high altitudes when convection and radiation become comparable 
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Heat Transfer Coefficient, C,, 


Ames Research Center 


. f222R?q(0)sino a0 ; ; 
a een’ ouA * C,, based on hemisphere computations 
es 2 ° Will be slightly different from full sphere (wake) 
E Cn anr30000 - Peaks at stratopause (roughly) 
J | . 
L(t os oe * C,, decrease in stratosphere due to 
s na exponentially increasing atmospheric density 
a — ¢ Discrete data curve fit in altitude (Z), velocity 
ue <a (u,,), and radius (R,,) 
Rectal sedinaio hs 
30 m diameter, Velocity variation 20 km/s Velocity, Diameter variation 
#V=12km/s @V=16km/s ®V=20km/s +D=1im *D=3m 4D=10m *D=20m BD=30m &®D=100m 
Ni Ni 
ae ts — 
7 3 3 a 2. 2 
fo) £ a ° £ 
Wipe = = Be ee eB S 3 om C= | Wi a Senne OO 2 38 
. / : 7 oe oa = is 
a s P : 2 8 g A f 8 g 
6 ° EY ae aft? 
° om S a S| S 
Meso. rie Strat ® FF o = __Meso. = Strat. __* a o = 
100 90 80 70 60 50 40 30 20 10 O as 100 90 80 70 60 50 40 30 20 10 0O 7 
Altitude, Z/km Altitude, Z/km 
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~ | Body+: Wake (L = 6D) 
Gas cap 
| I 


Temperature/K 


7 0 12000 24000 


50 ; T T 100 T T T T 150 ‘ T 1 
Axial coordinate/m 


30 m diameter, Velocity variation 


Luminosity" 
(Methodology from Stardust Mission) 


¢ Modest number of computations for sphere 
° Diameters: 1, 3, 10, 20, and 30 m 


¢ Wavelength range: 85 nm to 4 um 


¢ Radiance integrated over projected area 


20 km/s, Diameter variation 
*D=1m *D=3m 4D=10m *D=20m S8D=30m _ 


*V=12 km/s @V=16 km/s ®V=20km/s 
e 
B 
e 
a 
¢ 
o 
- oO 
o fol : 
Oo a 
Meso. Strat Zz 
c emg: mee eee | T T T T eo! heme T y T T "T 7 rt. r T T 

100 80 60 40 20 0 

Altitude/km 
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° Velocities: 12, 16, and 20 km/s 
¢ Stagnation pressures: 1, 10, and 100 bar 
¢ Lines of sight divided into 3 groups — 


nosecap, body, and wake 


1.E+06 1.£E+08 1.E+10 1.E+12 1.£&+14 


Emitted Power/W.sr* 


a 
e 
B A 
+ 
4 
Meso. Strat. 
a a an a a a a a a a ee a a 
100 80 60 40 
Altitude/km 


o8 


20 


Trop. 


1.E+06 1.E+08 1.£+10 1.£4+12 1.£+14 
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Emitted Power/W.sr? 


Multiple Bodies & Interactions 


¢ Currently no model/mechanism for fragmentation — various hypotheses 

¢ Supplement current knowledge with computations for idealized shapes 
e — ¢ 3D computations for various 
wea idealized shapes and arrangements 


) = in 
° Extraction of wake luminous energy 
° Extraction of aerodynamic/ i 
aerothermodynamic interaction forces/ :_ 


energies from computations 


dinate/m 
dinate/m 


. ow al . 
=; . Significant resources required!! raat fi 
0 20 40 “ a 3 20 40 
Axial coordinate/m Axial coordinate/m 
| Tale in i | le . 

E L b= £ 

2 2 2 2 

a a a a 

£ & £ £ 

z z z z 

So So So Se 

o* oe c* 6° 

uu uu uv uu 

® z z ® 

& 2 2 2 

tc t t t 

& & . 

> 2 o 

20 40 a 0 20 40 69 0 20 40 4 20 40 
Axial coordinate/m Axial coordinate/m Axial coordinate/m Axial coordinate/m 


All results shown here are for a velocity of 20 km/s and 30 bar of stag. pressure 
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ise Irregular Shapes 


* Conventional meteor physics assumes a spherical shape at entry 

¢ Computations on scaled versions of Asteroid ltokawa & NEO 2008 TC, 
¢ Itokawa is a dumbbell shape with “weakness” at the neck 
* NEO 2008 TC, most likely oriented in flight 

¢ Irregular shapes will require full (DOF analysis 
° Will require mass moments of inertia of the object — difficult! 


2x-scale NEO 2008 TC, 1/38-scale Itokawa 


ine 


9 (kWiern’) 
™) 45000 
37500 
30000 
22500 
15000 
7500 


2008 TC; geometry — Scheirich et al., MAPS, 2010 Failure likely 


7/7/2015 Page 14 


Lf Circling Back 
‘Ames Research Center 
¢ Can some of the modern computational analysis tools, used in design of 
heatshields of capsules, be used (repurposed?) for simulation of meteoroid 
entries? 
¢ Current limitation on entry velocity 20 km/s will be removed by improvement of 
thermodynamic and transport properties 
— Include doubly- and triply-ionized species (N2*+, 07+, N’+, and O%+) up to 50,000 K 
— Paper by Jaffe et a/. to be presented at AIAA SciTech 2016 in San Diego 
e Thermal response models (ablation/recession) being developed for silicates 
— Paper by Y-K Chen to be presented at AIAA SciTech 2016 in San Diego 
e Expansion of spectroscopic databases for silicates and metals under way 
— Effort led by C. Bauschlicher and A. Brandis 
e Radiative heating computations are not very efficient — need to replace process 
— Opacities for high-temperature air & stony meteoritic vapor 
— Account for radiation blockage by meteoritic vapor 
e Tighter coupling of analysis tools is required 
— Rapid recession is a hurdle, but not an insurmountable one 
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Supplemental Questions 


‘Ames Research Conler 
¢ How could the simulation tools and processes be enhanced/improved? 
¢ Development of material thermal response models, including multiphase flow 
e Physics-based models of fracture based on observation of recovered meteorites, 
fracture mechanics tools 
— Guidance from Ames Chief Engineer 
¢ How could these tools be verified/validated? 
¢ Simulation of well-known bolides (Chelyabinsk, 2008 TCs, ...) 
— Look to the meteor physics community to define test cases 
e Pathfinder experiments (see posters of Entry Technology Division) 
— ARC ballistic range or DLR wind tunnel for fragmentation 


— Arc jet for material response to convective heating and spectra of shock-heated air and 
ablation products from meteorites (or surrogates) 


— High pressure arc jet testing possible at AEDC 
— Shock tube for thermochemistry of shock-heated gases and “end wall” testing of 
meteorites for fragmentation 
¢ Airborne observation campaign(s) (see poster by Jay Grinstead) 
— To obtain flight data on fragmentation and bolide spectra 
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Final Thoughts/Opinions 


¢ How valid and useful is the assumption that a large asteroid is a sphere? 


e Entry vehicle heatshields are no longer designed by stagnation point environments 
alone — size and shape matter 


¢ How valid and useful are the “textbook equations” of meteor physics? 
e Meaningful for a single body, but perhaps not for a collection of objects, especially 
if there are interactions and a large loss of mass during atmospheric flight 
¢ How valid and meaningful is the mass loss equation of meteor physics? 
eA surface mass and energy balance seems to be the right way to go 
e Reformulation is necessary if multiple phases are involved 


¢ How do material properties scale from the lab to an exo-atmospheric body? 
e How valid is the assumption of isotropy? 
e Will require a new approach to structural analysis of porous media with internal 
cracks 
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“Ames Research Center 


Backup 
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Importance of Shape 


¢ For the same entry velocity and 
flight path angle, two bodies of 
different shapes and sizes, but 


Ballistic coefficient 


ASPECT RATIO = 0.8 ASPECT RATIO = 0.6 


(11.60 m x 9.28 m) (14.06 m x 8.43 m) identical ballistic coefficients, will 
eond/W.cm fly the same trajectory 
be ¢ Shape and size determine 
~ heating, and hence, mass loss 
72 ° Insufficient to use albedo inferred 
0 dimension to estimate the mass 
C,, = 0.938, A = 338.42 m?* C, = 1.013, A = 372.48 m? fora spherical shape 
ASPECT RATIO = 0.2 ASPECT RATIO = 0.4 e Itis important to know both entry 
(29.24 m x 5.85 m) (18.42 m x 7.37 m) 


mass and cross sectional area 

¢ Drag coefficient can be estimated 
easily using simple Newtonian 
impact theory 


¢ Videos of Chelyabinsk dust trail 


C,, = 1.192, A = 537.20 m? C, = 1.101, A = 426.38 m’ suggest a prolate ellipsoid shape 
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Ames Reseorch Cenler 


Entry 
interface 
(100 km) 
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M,, 

e= C, A, 
ae —u, siny 
dt " 
dy == cosy- 
dt U,, 

du, __ | PU, 

dt 2 £B, 
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Ames Research Center 


7/7/2015 


Energy 
Deposition (& 
Fragmentation) 


Simulation Process 


Entry State 
(Shape, Size, Spin, 
Material, Fracture) 


i i yc ey 


ENGINEERING-FIDELITY 


F 
al 
Wu 
= 
a3 
as 
© 
= 


—, Analysis Tools, I:" 
Sealeoeik tu Flight Mechanics 


¢ Single body; point mass (8DoF) 


e Runge-Kutta time-integration of equations of motion 


¢ Requires: 
Shape and dimensions of entry object Entry flight path angle (inertial) 
Entry mass Lat./Long. at atm. pierce point 


Entry velocity (inertial) 


e Modification to include mass loss 


e Time-varying heat transfer coefficient based on flow computations for 
spheres 


¢ Open Issues 
¢ Multiple bodies (due to fragmentation) and their interactions 
¢ Irregular geometries/shapes 
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J Analysis Tools, Il:" 
fees Roch Cone Aero/Aerothermodynamics 


e Flow (11-species — N,,0,,NO,NO*,N,*,O,*,N,O,N*,O*,e" — air model) 
e Eulerian frame of reference — flow past a fixed body (shape, size, orientation) 
° 3D Navier-Stokes equations 
— Turbulent flows use simple eddy viscosity model 
— Thermal and chemical non-equilibrium (rate processes) 
— Bow shock captured as part of solution and grid tailored to bow shock 
— Variety of surface BCs — radiative eq., recombination rate chemistry, ... 
¢ Radiation 
e Line-by-line spectral simulation; continuum included 
— Atomic line and diatomic spectral databases 
¢ Tangent slab model for radiation transport 


¢ Open issues 
e Multi-stage ionization required for velocity > 20 km/s 
¢ Coupling between material response and flow/thermal radiation fields 
— Solution turnaround time and computer resource requirements 


¢ Tangent slab transport model is somewhat restrictive 
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an Analysis Tools, IIl:" 
hans Novorh Cone Material Thermal Response 


e Multi-dimensional analysis of ablation 
e Detailed surface energy and mass balance 
— Ability to handle ablation of silica-based or pyrolysis of organic (phenolic) materials 
e Mass transfer handled through non-dimensional parameter called B-prime 
— B-prime tables (equilibrium) computed over a range of pressures and temperatures 
— Provide blowing rate (& composition) of pyrolysis products 
e Ability to handle recession of wetted surface 


e Requires (especially at temperatures above 298 kK): 


Mass density (and porosity) Elemental composition of material 
Heats of formation of constituent species Specific heats of constituent species 
Thermal conductivity Coefficient of thermal expansion 


¢ Open issues 
«Compositions of meteoritic materials vary 
¢eComplexity of thermochemistry of ablation products and shock-heated gas 
¢ Tight coupling between flow/radiation and thermal response is a challenge 


— This is required to assess radiation blockage by meteoritic vapor 
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—SJ Analysis Tools, IV:" 
foes Beach Cou Material Structural Response 


¢ 3D finite-element approach to static and dynamic loading; nonlinear 
analysis 
e Variety of elements; 6-noded hexahedral elements preferred 
— Easier data transfer between flow and structural meshes 
¢ Structural, thermal, and thermal-structural analysis 
e Anisotropy in material properties 
e Crack and crack propagation; strain energy release 


¢ Requires (particularly at temperatures above 298 k): 


Material density Poisson’s ratio 
Material moduli Shear moduli 
Yield strengths (compression, tension) Porosity 


¢ Open issues 
¢ Relating meteorite internal structure (incl. cracks & voids) to actual object 


¢Can codes based on the continuum hypothesis be applied to objects with 
large numbers of cracks and pores (micro-scale)? 
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- 
Entry Capsules vs Meteors, I: asap 


4 Reser Come Easily the Hard Part 


Capsule (Earth entry) Meteoroid (Asteroid) 


Size Less than 6 m diameter Varies from microns to (kilo)meters 
Entry Mass | Usually small for robotic missions | Varies from few grams to millions of kg 
Entry Type | Controlled; Ballistic or Lifting | Uncontrolled; Ballistic assumed 

Entry Velocity _ | Between 7.5 and 13 km/s | Hyperbolic (11-30 km/s) 

Entry Angle No steeper than -13° Varies between -90° and 0° (skip out) 
Flight Dynamics | 3DoF or 6DoF 3DoF (6DoF?) 

Aerothermal Dominated by convective heating Dominated by radiative heating 


The lower limit on size is dictated by airburst/terrestrial impact risk 
The upper limit on size is dictated by computing resources & solution turnaround 
The upper limit on entry velocity is dictated by gas-phase thermochemistry 
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Bow shock — 


yj Flow separation 
Boeniany aX \ | y— Shear layer 
layer ey 


;— Transition to turbulent 


Viscous —_ \ a, 
interaction \ : y~ Impingement (reattachment) 


Surface 
recombination —~_ 


Radiation —— 


ata 


Bit-t-telett lilelac 
ionization ’ 
(thermochemical S Reaction 
lelamueHillisselt ud] wt - retelalige) | 
Ss ; plumes 


-—_ 


Transition to — a 
turbulent , 
Ablation —“ ~ 
“~~ Control 
Shock-shock Interaction “ surface 


Figure courtesy R. Wheless & T. Horvath, NASA Langley Research Center 


Flight Mechanics 
¢ Ballistic vs Lifting entries 
¢ S3DOoF vs 6DOF 


Aero/Aerothermodynamics 
¢ Turbulent convective heating 
¢ Thermal radiative heating 


Material Thermal Response 
¢ Ablation & Recession 
¢ Melt vs Vaporization 


Structural Response 
¢ Pressure & shear loads 
¢« Dynamic loads 


High-fidelity tools are used to anchor faster engineering-fidelity tools via scaling laws 
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| Energy Deposition — Single Body 
Ablation 
coefficient 
Cy 
Rowe 


1, 1, 
I, =T|—Ou,, +1|—p,u,,Cp A O 
m 9) m 5 Patt D,, Am 


¢ tis the efficiency of converting kinetic energy (m,,u,,°/2) into luminosity 
—Should 4,, be the wetted area instead of X-sectional area? 

e Luminous efficiency is either specified or varied to match observed data 

e The heat of ablation, QO, a big source of uncertainty, shows up here too 
°Q is assumed same (8.08 MJ/kg) for stony and iron meteors, and everything in 

between 
° If C,, and /,, can be computed for a range of R,, (hence 4,,), p, (from p,,,.), and 
u,, of sohere, t can be estimated 


I 


m 


U, +] oC, An, 
2C, O 2 : 
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Ames Reseorch Center 


Stagnation line temperature/K 


600 


Heat flux/kW.cm? 
400 
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Shock-layer temperature 


P12km/s “~14km/s 


*16km/s —-i8km/s 20 km/s 


0 50 100 150 200 250 300 350 
Stagnation pressure/bar 


Surface heat flux distribution 


(u,, = 20 km/s, p,,.. = 30 bar, R,, = 15 m) 


—Convective —-Radiative —Total 


0 10 20 30 40 50 60 70 80 


Angle from stag. point/deg 


Flow Characteristics 


Free electron mole fraction 


@12km/s “14km/s ©16km/s —18km/s -®20 km/s 


03 O04 O35 


Electron mole fraction 


0.0 O1 0.2 


0 50 100 is0 «6200 0 «=6©250 063000 350 
Stagnation pressure/bar 


e High stag. pressure => low altitude 

e At 20 km/s nearly fully ionized flow 

¢ Boundary-layer is fully turbulent 

e Surface imperfections may enhance 
convective heating locally 

¢ Surface heating completely dominated 
by shock-layer radiation, especially at 
large length scales 
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Ames Research Center 
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30 m diameter, Velocity variation 


(185, A:0)da)a0 


Luminosity: Gas Cap 


¢ Gas cap luminosity is simply total energy 
radiated into a hemisphere enclosing the 
bow shock 
¢ Wavelength range: 85 nm to 4 um 


@V=12km/s ®V=16 km/s ®V=20km/s 
B [ 
° es = 
a ¢ [ 
° his F 
° a b 
o 
i 
Oo} 
Meso. Strat. = 
i T T | ee meme | T a T T : ee 1 + io ot co 
100 80 60 40 20 0 
Altitude/km 
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Emitted Energy (P°*,, ,,) /W.sr* 
1.E+06 1.£€+08 1.£+10 1.£+12 1.£+14 


* Can be tailored for any passband filter 
¢ Conversion of energy to magnitude over a 
flight trajectory provides a synthetic light curve 


20 km/s, Diameter variation 


+D=1m *D=3m 4D=10m *D=20m ®D=30m ®D=100m 


e e¢ ® 
es 
D Q = 
a--3--% 
5 ra ° ke ck 
a = . ‘ 
° 4 
a 
co A 
a 
oO 
A 
Meso. + Strat. 
a. ik. t i 
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TMM} 
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Circling Back, II 


‘Ames Research Conter 


¢ Can we build or develop, across various classes of meteoroids, models for: 
e Material thermal response — ablation (vaporization vs melting) 
— Attempt at estimating enthalpy of ablation 
— Intense CO, laser heating (5 — 20 kW/cm?) of meteoritic materials (Tamdakht, ...) 
— See poster by Stern ef a/. in this workshop 


¢ Material structural response — fragmentation (internal structure) 
— Conversion of cracks/fractures in meteorites (falls) to unit problems for structural analysis 
— Design of Experiments approach to “unit problem” of matrix+inclusions 
— Influence of structural properties, including material anisotropy, of meteoritic material 
— See poster by Agrawal et al. in this workshop 


e Energy deposition along the meteor trajectory in the atmosphere — light curves 


— Attempt to reconstruct Chelyabinsk trajectory using constant C,, and time-varying C,, 
derived from flow and radiation computations 


— Fragmentation (Borovicka et al., Nature, 2010) imposed at 40 and 30 km altitude 
— See poster by Allen ef al. in this workshop 
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